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ABSTRACT 

The use of low thrust devices with continuous and discontinu- 
ous thrust programs is investigated to determine whether or not a 
discontinuous thrust program will provide a greater payload and 
structure mass fraction, where the program must impart a specified 
amount of energy to the satellite in a specified total time, starting 
from a circular parking orbit. A discontinuous thrust program is 
developed, based on the elliptic orbit, using a perturbation analysis, 
and a series solution is obtained which permits investigation of this 
program up to angles of one radian either side of the perigee of the 
elliptic orbit. These results are compared with a continuous thrust 
program which gave a spiral orbit. Under certain conditions, where 
storage batteries must be carried as part of the payload but are 
available for use during the thrust program, the discontinuous thrust 
program is found to provide a greater payload and structure mass 
fraction. Further investigation at angles greater than one radian 
either side of perigee appears to be warranted. 
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I. INTRODUC TION AND SUMMARY 



Once in a parking orbit, low thrust devices may be used to 
alter the path of a satellite.. In particular, it is anticipated that one 
use of the low thrust device will be to add energy to the satellite for 
an escape or a near escape condition. The exact path and thrust 
program must be determined by the mission requirement, but in any 
event, the system should be optimized for efficient launching. 

The problem of optimizing an arbitrary mission is a for- 
midable task to say the least. To simplify matters, it is assumed 
that the mission is to impart a specified energy in a specified total 
time to the maximum fraction of payload and structure mass. To do 
this, one must make efficient use of the power supply and propellant. 

The usual variation procedures may be used to derive the 
differential! equations for optimum direction and magnitude of thrust 
when continuous thrust is employed, as shown by Irving in reference 
(1). He has shown that a constant acceleration is optimum for the 
gravity-free case. No analytic solution has been found for the dif- 
ferential equations in the central force field, but Irving has made a 
numerical study and found the acceleration to be practically constant 
for various escape programs. 

In reference (2), Casey has shown the distinct advantage of 
energy addition by thrusting near the perigee of an elliptic orbit in 
such a manner that the perigee remains constant. His thrust pro- 
gram is formulated to require that the perigee distance remain con- 
stant throughout the thrust cycle, and results are obtained by 
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numerical procedures. The reason for thrusting near the perigee is 
to take advantage of the greater velocity, since in any elliptic orbit 
the velocity is greatest at the perigee and since the energy addition 
rate is proportional to the velocity. Further, if the perigee distance 
is kept constant, then the perigee velocity will increase as the eccen- 
tricity of the elliptic orbit increases. 

As pointed out by Casey, there will be a period of coasting 
between each thrust cycle which will make a considerable contribu- 
tion to the total time. For this reason, we cannot say simply by in- 
spection whether the continuous or discontinuous method of energy 
addition will provide a greater fraction of payload and structure 
mass, where both methods are required to impart a specified energy 
in a specified toted time. 

To investigate this problem, we shall use the idea of re- 
stricting the perigee distance, but only at the beginning and end of 
each thrust cycle. This permits added freedom in an effort to gain 
more energy. Since we are only considering low thrust devices, a 
perturbation technique is used with a simple variation to develop the 
first perturbation integral equations for the maximum change of en- 
ergy during one thrust cycle and for the Lagrange multiplier neces- 
sary to satisfy the above perigee restriction. The equations apply 
for any period of thrusting up to one complete circuit of 2tt radians, 
but could not be integrated in closed form. By making a series ex- 
pansion of the integrands, the equations may be integrated term by 
term to obtain results useful to an angle of one radian either side of 
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the perigee. While this does not permit a complete comparison of 
these two methods of energy addition, it is possible to draw definite 
conclusions for angles less than one radian. 

After obtaining equations for one cycle of energy addition, 
the effect of a large number of cycles is obtained by integrating the 
changes of the various elliptic orbit parameters. 

To compare these results, an approximate solution for the 
spiral orbit, using alow thrust device, is obtained* Then a compari- 
son is made of the payload and structure mass fraction, where both 
methods are required to fulfill the mission of imparting a specified 
energy in a specified total time. 

Finally, since the discontinuous thrust program could 
generate energy while coasting, the use of storage batteries to aug- 
ment the power available during the thrust cycle will be considered. 

In one case, the batteries will be considered as part of the power 
supply; and in another, they will be taken as part of the payload, but 
available for use during the thrusting program. 

When thrust is only permitted up to an angle of one radian 
either side of perigee, the discontinuous thrust program was found 
to have no payload and structure mass fraction advantage unless bat- 
teries are employed. When the silver-zinc battery is considered as 
part of the power supply, there is an area of marginal advantage due 
to the low specific mass of this battery, but its poor recharging re- 
liability rules out its present use in such a system. 

The nickel-cadmium storage battery, which is very relia- 
ble, provided an area of definite advantage for the discontinuous 
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thrust program when the battery is considered as part of the payload. 
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II. PRELIMINARY ANALYSIS 

We shall first write the equations of motion for the satellite, 
making the usual simplifying assumptions of a very large, central, 
spherical mass and take the coordinate system at the center. We 
shall further assume motion in a plane, keeping the thrust in the 
plane of the orbit, and neglect the effects of other bodies. Then, 
using plane polar coordinates , as in fig. 1, and the dot notation for 



time derivative, we have: 

n„ =7 n ( 2 . i) 

Tty* = / fl +■ j ft & 

7f„ = r f/i -ne z ) + 7(7fiU n *®?) < 2 - 3 > 

Equating equation 2. 3 to the acceleration due to thrust, Ck , and grav- 
ity, we have: 

Tin-ne 1 )^ '«(?»() - 

+ j (o. S»N(|34-cJ))J (2. 4) 



where R 0 is the radius of the central mass, and $ 0 is the acceler- 
ation of gravity at fi = R e . Since a /^° is small for a low thrust 
device, we shall first take 0.-0 and obtain the usual solution in 
dimensionless form. The thrust will be considered later as a per- 
turbation. Choosing O as the independent variable, and letting ( )' 
be differentiation with respect to & , then with U as the dependent 
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variable and fx as the dimensionless angular momentum, we have: 

M R " 

u ~ 7 f (2.5) 



n z e = ]fi7^ Ro ft 

but when a = O , j- ( fl Z & ) = 0 

^U 0 • 

© - ]/^7 fi 0 u 2, 
ri = - IjfrRc -ft, u' 
n = - $* u 



( 2 . 6 ) 

; therefore & = constant = 

(2.7) 

( 2 . 8 ) 

(2.9) 




(2. 10) 



Substituting into equation 2. 4, taking O. - O , and simplifying, 
we have: 

U" + U = Jz ( 2 . 11 ) 

A solution may be written 

u = = 

where C 0 is the eccentricity and ©/> is the angle of the perigee 
position* The non-dimensional time. A' T* , for any portion of the 
orbit is given by: 



i + e o cos (e - © P )j 



( 2 . 12 ) 
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de 



(i t e 0 cos &y 



B, 



(2. 13) 



using equations 2. 7 and 2. 12. The integral can be evaluated to give: 



AT= 



(i 




t 

I +e o 2; 



(l + e.) -h(i-£ B )t4d z £ 




In equations 2. 14 through 2. 26, we shall essentially follow the 
work of Irving, reference (1), since his approach provides a conveni- 
ent parameter for computing the mass fraction of the satellite. 

The mass of the power supply, , is assumed to be pro- 

portional to the power supplied, P : 

= KP (2.14) 



where K is the specific mass of the power supply* It is further 
assumed that the maximum power is utilized when thrusting, appear- 
ing as kinetic energy in the exhaust: 

P = j(- *”) c 2 ( 2 . 15) 



where C is the exhaust velocity, vvi the mass of the vehicle at 
time t , and -vv? is the propellant flow rate. We also have the 
acceleration of the satellite due to reaction: 

thrust YV1 C 



a = 



mas s 



Vn 



( 2 . 16 ) 



Using equations 2. 15 and 2. 16, we eliminate C : 
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(Z • a. _ -vvi 2 ' [ C • C ] / /yy\ z 
2P - Yn C Z 



y*2_ _ / i \ 



(2. 17) 



Equation 2. 17 may be integrated to obtain an expression for mass as 
a function of time: 



_L = J_ + -L (*-£1 Jt 

M M c 2 J p 



(2. 18) 



where is the initial mass of the satellite. If we let be 

the propellant mass, then the "burn out" mass. Mly , at time 
is given by: 

M b =M 0 -M, (2 ‘ 19) 

Substituting in equation 2. 18 and using 2. 14, we have: 

^ = I ^ JL n. (\‘jt < 2 - 20 > 

M k 2 m„ ) 

0 

To optimize the power supply, we write the relation between 
the payload mass, M, , the structure mass, M$ , and the other 
masses: 

M, «■ M, + = M„ (2.21) 



Further, let 



y 2 = ^ 



(2.22) 



Then, manipulating equations 2. 19, 2. 20, 2. 21, and 2. 22, we have: 



M, +M g _ M w f 
Ho 



I 



M 0 [(Mw/Mj + y 



(2. 23) 
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There is an infinite number of programs for £6 which will keep Y 2 
constant while fulfilling the mission as is varied. Thus, 

taking the derivative of equation 2. 23 equal to zero, we have the 
maximum of ( M , fl. with respect to when: 

= Y - Y Z (2-24) 

Substituting equation 2. 24 into 2. 23, we have: 



M, + M s 
Mo 



(i- rf 



(2. 25) 



Using equations 2. 21, 2. 24, and 2. 25, we have: 

= Y (2. 26) 

ML 

From equation 2. 26 we see the range of Y must be zero to one, 
since Mp/M© must lie in this range. Further, for a particular 
mission, we see from equation 2, 25 that the structure and payload 
mass fraction will be greatest when Y is a minimum. Alter de- 
termining the minimum value for Y , the optimum power supply 
mass is given by equation 2, 24. 

We are now ready to investigate equation 2. 22. In reference 
(1), Irving has shown analytically that the acceleration must vary 
linearly with time in gravity-free space, if Y is to be a minimum. 
Further, when the mission of specified energy in a specified time is 
considered, the optimum acceleration becomes constant. For the 
case of a central force field, Irving made a numerical study using 
continuous thrust, and found the optimum acceleration to be practi- 
cally constant. Thus, the spiral orbit with constant acceleration is a 
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n near optimum 11 solution which will provide a fair comparison be- 
tween continuous and discontinuous thrust programming. 

If we consider a discontinuous thrust program applied to an 
elliptic orbit, it seems clear that the thrust should be applied near 
the perigee, where the velocity is greatest, to impart the greatest 
energy® If the thrust is parallel to the velocity vector, we will im- 
part the maximum energy during any given thrust cycle. Such a pro- 
gram will allow the perigee distance to increase, and with this in- 
crease there is a decrease of perigee velocity, thus decreasing the 
energy imparted during a subsequent given thrust period. To obtain 
the advantage of thrusting near the perigee, but still keeping the 
same distance for a velocity advantage in subsequent cycles, a pro- 
gram must be found that imparts maximum energy during any given 
thrust cycle while providing no net change in the perigee distance. 
After obtaining the restricted perigee equations for one cycle, the 
programming of thrust time will be considered and Y will be com- 
puted for a large number of cycles. 

After developing the restricted perigee equations, we shall 
obtain an approximate solution for the spiral orbit with continuous 
thrust to compare with the restricted perigee thrust program* Fi- 
nally, we shall consider the use of batteries with the restricted peri- 
gee thrust program. The batteries can be charged while coasting to 
provide added power when thrusting. In one case, the batteries will 
be considered as part of the power supply, and in another case they 
will be considered as part of the payload. 



- 11 - 



III. PERTURBATION EQUATIONS 



When Cl is not zero, we have a small non-dimensional ac- 
celeration: 



a 






(3. 1) 



since we are considering only low thrust devices. This suggests the 
use of a perturbation analysis for each cycle. 

Since $ is no longer constant, we must rewrite equation 

2. 9: 

fi - - U z (-k'u'-h& u") ( 3< 2) 

and for the j component of equation 2. 3 we also have: 



-L 

n 0 

By using fig. 1 and simple trigonometry, we write: 

-h 



(3.3) 



SIN(Cp+p) = ~==\ 



(3. 4) 



and 



p) = 



n'coflf) - US/AJtp 



(3. 5) 



Substituting equations 2. 5, 2. 10, 3. 1, 3. 2, 3. 3, 3. 4, and 3. 5 into 
equation 2. 4, we have the radial equation: 



U + u +- 



fV 

-k 



= 1? + e 



U.' COS (f> f 

it z u z l)u* + u /a 



(3.6) 



and the tangential equation 
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o / _ / -u'si»<j> + \ 

~ U 3 & j 



(3.7) 



We get some simplification in equation 3. 6 by substituting for ■£/ 
from equation 3. 7: 



u" + U = + £ 



y u z -h u' 2 



2 , .3 






7 A/ <p 



(3.8) 



For the perturbation solution, we take $ 0 constant and U 0 as 
given by equation 2. 12 to write the assumed solutions: 



and 



U - LA c 4- £ U , + £ + £ 7 U y 4- 



-fa — 4- £ , + £^& x 4- £*$ + 



(3.9) 



(3. 10) 

where U, - U^~ ■ * * = = ' ' • — O at O t , the angle where 

the thrust for a particular cycle is started. Substituting equations 
3. 9 and 3. 10 into equations 3. 7 and 3. 8, we obtain the first perturba- 
tion differential equations: 



j / __ -u'sisJ<p + U c c. os <f> 

' ’ ' 4 . ul 



(3. 11) 



and 



LA. 4- U . = 



- li 



, ifu/ 4- ul z7 „ M 

4- ■“ 1 — S/AJ Cp 

^ < u r 



(3.12) 



For a particular cycle, where the thrust is started at angle & t and 
G - 6, ~ Z IT , the solutions may be written: 



Ve) = ( 

e, 



(3. 13) 
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and 



© 



U/9) = j S iaJ(B-9) 

e, 






(3. 14) 



For a unit mass, the energy is given by 

E = kinetic energy + potential energy 

^Ro 

_ + pR'C I-U) 

= f.t.[ I - Ytf ] 



(3. 15) 



by using equations 2. 5, 2. 6, 2„ 7, 2. 8, and 2. 12, and where the po- 
tential energy is considered zero at n - . By substituting equa- 

tions 3„ 9 and 3„ 10 into equation 3„ 15, we obtain the first perturba- 
tion energy change: 



Al z_ _ 
■ 



- + (tfu o -l)U, + M,' 



(3. 16) 



The perigee distance ratio is obtained by evaluating equation 
2. 12 at 6 = e ? ; 

R 



' O 



+- £< 






(3. 17) 



We could require that this quantity remain constant during the thrust 
cycle, as done by Casey in reference (2), but we should have a gain 
in the energy imparted if we require that the perigee distance ratio 
be equal before and after the cycle. We shall pursue this latter 
method of perigee restriction and will demonstrate a slight gain of 
energy over the method of continuous perigee restriction. 
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The eccentricity may be evaluated before and after the thrust 



cycle: 



e 



'(urir + w 



-A 4 



(3. 18) 



by using equation 2. 12. Then, substituting equations 3. 9> 3. 10, and 
3. 18 into equation 3. 17 and equating the results before and after the 
thrust cycle, we obtain the first perturbation perigee restriction: 



pr = o = ^ - ?‘0 



-h 






(3. 19) 



evaluated at & 2 , where is the angle at the end of the thrust 



cycle. 
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IV. OPTIMIZATION OF THE RESTRICTED 
PERIGEE CYCLE 



To find the thrust angle which will impart maximum energy 
during the thrust cycle and still satisfy the perigee restriction given 
by equation 3„ 19, we shall take a simple variation, using the La- 
grange multiplier, equal to zero. We write the non-dimensional en- 
ergy, equation 3. 16, plus the perigee restriction constraint, equation 
3 0 19, with the Lagrange multiplier, X , where A is a constant to 
be determined, and, assuming 6 = constant during the cycle, take a 
variation on <p equal to zero: 



s 



4 E 









+ APR 



= 0 



(4. 1) 



For convenient notation, we rewrite equations 3* 16 and 3, 19 as 



~~~ = /U, + +■ Cuf (4 .2) 

t K o 

and 

PR = 0 = A% + Bu, + Cu,' (4 . 3) 



where A , A / , 13 , and C- are shown in equations 3, 16 and 3. 19. 
Here, we note, by talcing the difference of equations 4. 2 and 4. 3, that 

rf*r <*•*> 

is the energy change for the thrust cycle. Substituting equations 4. 2 
and 4, 3 into 4, 1, we have: 

(A + 1A')W, +(i + >)[6Su, + CSu/] = 0 (4.5) 

Before taking the variation, we must first integrate equation 3„ 14 by 
parts to obtain the more useful form: 
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© 9 

U,(e)4s«(S- + -| j\[coi(e-f-)-ij E'lWJ'f (4. 6) 

l “• i. 



Then, taking a derivative of equation 4* 6 with respect to © : 



a 



u/£>) 

d «„ Me 



& 

|}f 



(4.7) 



Using equation 3. 11, we note that: 



©fc 



^2 



s W =( ) sift rut = 



©, 






-w! Coi ■ <p - <,!AJ(f> 

t. u ° fufTu^ 1 



$<!>*? 



(4.8) 



Now, taking the variation of equations 4 0 6 and 4 0 7, using 4. 8, we 
substitute into equation 4 0 5 to obtain: 

S [«F«- +JPR ] ~S & i c u! * c ‘° >( O t-rjc-e*# ~ 



~7~" + ^z-“ ~ -y- S I/O fa- */*'<? + u' cos ^jj 



(4.9) 



For equation 4. 9 to be zero, we must require that the integrand be 



zero, so that 
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= - u'c°‘4>[ Aili' *^-(cpi(^-f)-i) - 44 - 



Z a Z 

Uo + Wp 






SrnJie^- V) + C c ^ r ($ 2 .-^) 



(4. 10) 



Substituting the expressions for A , A' , B , and C from 
equations 3. 16, 3. 19, 4. 2, and 4. 3, using equation 2„ 12, and taking 



Qp — O t equation 4. 10 becomes 



"ttXN d) = 

^ j + e 0 tose\_l + et + ze.& os e + zhe a (aos&- i\ 

where A is a constant to be determined. It will be shown that the 
shift of Of> is negligible for the symmetrical thrust cycle, so that 
a choice of zero for Qp presents no problem in going from one 
thrust cycle to the next. We now have the thrust direction program 
given by equation 4. 11. 

Substituting the expressions for A and A and equations 
2. 12, 3. 11, and 4. 11 into 4. 4, we have: 

A E 



(4. 11) 



£ R.% 



1 (i i-ej’i i + cw^f [ \ + e.} + z e. cos v + z ^e B (cos^-i)] - edn-e^ XeoSw f 






J ( 1 + e?+Ze*cos<r ■hZ’he„(cos + 

e ' 1/4 fefi l + e H* 

(4. 12) 



Also, substituting equations 3, 11, 4.6, 4.7, and 4. 11 into 3. 19, and 
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using 2. 12, we have: 

PR = o = 

_ r \e B (\j-e 0 ?snfi' ■hZ(He.cosH') l ll + el+ze.<ux t i' i-ne^co* * -i)\(cos¥-i ) ^ 

J ( l +■ e,c.os y) 3 fi + ef +ze, a>s v Ilf / ~+e 0 7^f[\ ~teT^ze c ~c°s V JITs V- 7) J 1 T ' 

II + retcn-ej Wr 



(4. 13) 



From equation 4 0 13 we determine the value of the constant A<5 0 
for any particular cycle, since always occurs with S 0 in the 

two equations 4* 12 and 4. 13. 

Since it is not possible to integrate the above equations exact- 
ly, we must use numerical methods to carry the investigation to 
large angles. This would be fairly simple for one cycle, but when 
one considers a large number of cycles, a computer program appears 
to be mandatory. We shall not pursue this route, but instead, we 
shall present an approximate evaluation which is useful to angles of 
one radian. 
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V. APPROXIMATE EVALUATION FOR ONE CYCLE 

To make an approximate evaluation of the optimum restricted 
perigee equations for one cycle, we shall make series expansions and 
integrate term by term. First, using equation 4. 13, we may rewrite 
4. 12 as: 

P ri[i + e*+ze.cciV)(H-e 0 cos<i‘fln-e^cie.c 0 sr+Z\eUcos</'-o] ^ 

_ _ Q Pfc +■ \ 

eMt J DENOMINATOR 

® 0. (5. 1) 



since PR - O by the appropriate choice of As o and where 
DENOMINATOR is the denominator of the integrand in both equations 
4. 12 and 4. 13. 

Now consider the function 

G> = ^ ^ ^ (5. 2) 

e, 

where 




")} 0 + e 0 coiH'f[(H-et+ze,cos i i')b 

( l + e 0 cosv ) 3 ]) i Tef + ~t e a co ~s v 



(5.3) 



We notice that 

_j_ ' S 6 _ _l_ ^ 

a ~ c„n 

b-i 




e, 



b = l 




JV 

y--i 



has an integrand identical with equation 4. 13, and 




(5. 4) 



(5. 5) 
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has an integrand identical with the non-vanishing term of equation 
5. 1, since PR - 0 . Thus, if we perform the integration of equa- 

tion 5. 2 , we may evaluate /\ e 0 and the increase in energy by 
taking the derivatives indicated in equations 5. 4 and 5. 5. 

Expanding in series, we collect the coefficients of the various 
terms of $ and integrate. Noting that the result is symmetrical 
about 0=0 , we write 

J_r _ f/q + N> 0* , Na ef Ns g 7 ? 

2^ (l+e m )l z 0+e o f G "(l +e 0 fi zo io ,080 * ' ' ‘j (5.6) 

where 

N, = aef +e 0 - zr + a - 2 ( 5 . 7 ) 



= I6<2„ 4 + I9ei-e„ -e. -zz^el -22 <re 0 + z <r + 20 T 2 ef 4 

+-J4-T z e o - 4 I 2 <r 3 - 3 T 4 (5.8) 

N 3 = + I03eef +z4ze?-it>ze? - tee? +ze 0 -i+qref-WL re ?- 

-17 re* -t-3 +4 re 0 ~4r -h izrz r z e* + 3 & 7 $<r z e<? -*■ 28 oz r z e z - 
-Sza. (r z e 0 4 - 72 <r 2 4 u*t> r 3 e? +Z& 4 o r^e B -300 r? - 4 zo r^e* - 
- 93of%, 4 24 or*- 540 -h do (T 

(5. 9) 



and for convenience we have used 

r = 

b (5. 10) 

Now, taking the derivative of equation 5. 6 as indicated in equation 5. 4 
equal to zero, we satisfy the perigee restriction of equation 3. 19 and 



solve for A e o : 



- 21 - 



^\-rA- 



wse?- st.'tzel ~jo&s^e^ -d33&e 0 - 3^92 &1 

(i+ej 4 i~g^&oo 



(5. 11) 



Here, we note that equation 5. 1 1 would be much more complicated if 
a non- symmetrical thrust program were used. However, since we 
can see from the form of equation 5. 6 that this will have no energy 
advantage, a symmetrical thrust program is chosen to simplify the 
results as much as possible,, 

Similarly, taking the derivative of equation 5. 6 as indicated in 
equation 5. 5, and using equation 5* 11, we have the energy for one 
half of the cycle: 

i =i ±1 f s _ a-ze.) el + (7 - ue^ + ize?) ef 
Z £$° R ° \+e 0 \ 2 - (i+e.) 6 (l+e 0 ) z izo 

f 7 

zizo e„ + 5i9o ef+2ez4-e?-Z504-el - 15 -os e^~ii34-e 0 -ni I e z 1 
TT+elf FoToo ) 



(5. 12) 



To determine the advantage of this method of perigee restriction 
over that used by Casey, the energy equation given in reference (2) 
was expanded in series and integrated to the same number of terms 
as we have in equation 5. 12. The two energies are identical to the 
last term, but if we compare the last term of the constant perigee 
energy expression: 

__ 2720 el -hsiooef+u zSet -zoooel - zssoel - 4Z30g 0 - 2,335 el 

t ^ ’ ~~ " (5. 1 

0+e,) 6 5o4oo 

with the last term of equation 5. 12, we see that for values of S c 
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between zero and one, the constant perigee restriction yields less 
energy. 



To find the changes of the elliptic orbit parameters U > 'h , 
and <2 , we note from equation 4. 4 that 




l_ 21 £ 

(A - A / ) 6 ft" 



(5. 14) 



where the energy is given by equation 5. 12. Substituting the expres- 
sions for A and A / and neglecting the higher perturbations since 
€ is small, we have 



A & — £ 



-ft! 

fiTij 2 <f 0 R o 



(5. 15) 



Now, from the perigee restriction, 
and after the thrust cycle: 

I + _ I f g, 

t. 2 “ 



we may write equation 3. 17 before 



(5. 16) 



from which we solve for the change in eccentricity: 



A e = 2. 



0 +e.) 



A E 

$ oRa 



(5. 17) 



using equation 5* 15. To see if there is a shift of perigee, we assume 
that the equation for U , at the end of the thrust cycle, has the form 
of 2. 12 : 



[i +(e o ^e)co^(& z -A0 p j] 

U= (t.+Ai) Z 

where we have provided for a shift of the perigee with the term 
Then, using equations 2. 12, 3. 9, and 5. 18, we can write 



(5. 18) 
A &p * 
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(t. + 



(5. 19) 



and 



U 



/ 



[ 




(5. 20) 



( lo 



Substituting these two equations into 4. 3, with the expressions for 



( Cos (G^-A^p) -cose^cos 6 Z -o.b p )-s/aj - O + 0(€) (5.21) 

Thus A&f> is at most of order £ . Using this result, we may 

evaluate the terms of order £ in equation 5. 21, using 5. 15 and 
5„ 17, and we find that they are zero. Thus, there is no shift of the 
perigee within the magnitude of the first perturbation as would be ex- 
pected with a symmetrical thrust program. 

Before proceeding with the discontinuous thrust program, we 
shall obtain the approximate solution for the spiral orbit with thrust 
tangent to the velocity vector. For this case, is zero, and equa- 

tion 3. 8 become s 



j\' , 0 , and £ , we obtain the result 



n u" + U 

From equation 3. 7 we have 



(5. 22) 




(5. 23) 



Talcing twice of equation 5. 23 equal to the derivative of 5. 22 and ex- 
panding, we obtain 
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where we have assumed starting from a circular orbit with 
and (u'/U) much less than one. Integrating from O to 
have: 



(5. 24) 
( u"/la) 

& we 



u 2 = ul- tee -z 



HU -h 4 £ 



U. 

u 



~ € (u) Z+ ‘") el6 (5.25) 



To determine approximately the small quantities not yet integrated, 
let us take for the first approximate solution 



ul = ul -4 £0 

which is merely the first result from equation 5. 25. Then, 

derivatives of equation 5. 26, we have approximately: 

U. ' _ Z_6_ 

u Ui 

j/' 4 6 2 

U " Ul 



(5. 26) 
taking 



(5. 27) 



(5. 28) 



UM»= - 

Ua 



(5. 29) 



Substituting equations 5. 27, 5. 28, and 5. 29 into 5. 25, we obtain the 
result: 



— ul + ZZ 4^ U/ 2, U/ 2 ~j (5» 50) 

where is given by equation 5. 26. We may now use the second 

term of equation 5. 30 to choose values of 4 and U , assuming 
U 0 is given, to insure that our results lie within a certain accept- 
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able error. We should also note that ( 6/ U A ) must be much less 
than one if equation 5. 30 is to have meaning. 

To determine the time, we use part of equation 2. 13 and sub- 
stitute equations 5. 22 and 5. 30. Then, using equations 5. 26 and 5. 28, 
we obtain 



A 777 & 

? fjU+U , r , ..ill 

( i + uh(«) 

0 0 



' ") ede 



6 * = fz? - iK ~ 






37 



u 



Vz. 



22 

u* u* 



4-0 "1 

uH 



(5.31) 



Here, we see that for specific values of and G , we will have 

a larger error in equation 5. 30 than in 5. 31, since is less than 

IA 0 . and U 0 is less than one. 

To determine the energy input, we substitute equations 5. 22, 

5. 27, 5. 28, and 5. 30 into 3. 15, and take the difference between the 
final and initial values: 



A £ 



l 



■ Rc 



- £ ( Uo 6 ( Ul 6/ 0 2 ^J] 



(5.32) 



where again the error may be estimated by the second term. Of 
these three equations, 5. 30, 5. 31, and 5. 32, equation 5. 30 will have 
the largest error as U&_ becomes small, and it should be used for 
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choosing values of 6 and U a . 

We have reached the point where we must decide how to pro- 
gram the cycles for the restricted perigee. To do this, we shall 
choose a ratio of (Ro/ft) — 0.9S to be the maximum value for U 
For the case of the earth, this corresponds to an altitude of 208 miles 
and is adequate to insure negligible atmospheric effects. 

If we form the quotient using equations 2. 13 and 5. 12, 



specific power = 17— - - 

( C K o 



A_E_ 

A'T 



then expanding A'T' in series and performing the division 

» = -L- £ + + ... ? 

( f,. I i +e ° e (i j 



(5.33) 



(5.34) 



This equation is plotted in fig. 2 for various values of <£ „ , with 

corresponding values of (E / R* ) for P 0 — 0.9b f where 
Po is defined as 



P. 



I + 



(5.35) 



Also plotted in fig. 2 is the specific power for the spiral orbit for the 
same values of (E / R &) . We can see that when the eccen- 

tricity is small, there is not much difference between the two methods 
near the perigee, but the discontinuous thrust program becomes inef- 
ficient if we proceed too far from the perigee. This is because the 
thrust must be directed away from the velocity vector to satisfy the 
perigee restriction. Thus, for a given energy increase per cycle 
when the eccentricity is small, £ must be larger with discontinuous 
thrust, since the thrust is not on for the full Z/lT radians. There- 
fore, 23 Y is larger, and the payload and structure to initial mass 
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ratio is reduced, as seen in equation 2* 25. If we manipulate equa- 
tions 2. 13 and 3. 15, we can obtain the time for 2 7T radians of trav- 
el in terms of the total energy: 



which holds for both elliptic and circular orbits. Now, if the mission 
is to be accomplished, both methods will have about the same number 
of cycles, where in the case of the continuous thrust method, cycle 
means £ IT radians of travel. Thus, we want to impart about the 
same amount of energy during each corresponding cycle. 

We can also see in fig. 2 that for larger values of eccentricity, 
say 0. 75, the specific power is much greater for the elliptic orbit* 

In this case, for a given energy increase per cycle, /\Y for the 
spiral orbit will be greater. 




(5.36) 
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VI. APPLICATION TO MULTIPLE CYCLES 

To find the total energy input, total time, and Y , we must 
make a summation of the results for each cycle. However, since we 
are assuming a large number of cycles, with small changes for each 
cycle, we may write integrals for the summation operation and ob- 
tain our results by integration. The thrusting time, A t , is taken 
as constant for all of the cycles, and the variable of integration is 

Je . 

We must first express the cutoff angle, 0, , in terms of <£ 

because it must be varied slightly to keep the thrusting time constant. 
Expanding equation 2. 13 in series, we have, for the symmetrical 
orbit: 






e 

0 + e) 



& z 



e0z_8e_) 

6o(i+e) z 



G z + " 



Then, if we let 






A ( 



( 6 . 1 ) 



( 6 . 2 ) 



we can write approximately 

where ji is the initial value of Q z . Care must be used in 

choosing 'fe so that ©^ is loss than one at fi = | \ otherwise, 

the series expansions will not be good approximations. 
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Substi tuting equation 6. 3 into 5. 17, using 5. 12, and replacing 
the finite "delta" quantities with differentials, we obtain 



lU<4l 



7 ae 3 - 29<? 2 -32£ -// .47 . 

+ -ft V de 



4 izl' £ ' f i +■ e 3 go (i ->-e) y/z 



(6.4) 



The term dt is in effect one, that is, one cycle, and (<$&■/ 6) is 
the small change of eccentricity during that one cycle. To perform 
the integration, we shall take £ as constant to obtain: 

( 6 . 5 ) 

2 - 
- L 

where <2/ is the initial eccentricity and <£?£• is the final. N is 
the number of thrust cycles, and the total thrusting time is given by: 

£T t = € Nat ( 6 . 6 ) 

We find the energy input very simply by replacing the finite "delta" 
quantities in equation 5. 17 with differentials, using 5.35, and inte- 
grating to obtain 






lezH-iszee-s?^ 2 -* z ice 7 0 ^} 
~Z7oo jJTe^ * 



E 



■P 





(6.7) 



Then, using equations Z a ZZ, 3, 1, and 6 a 6, 



r 



a _ 






' — /V 2J T 



( 6 . 8 ) 



To find the total travel time for 27T radians, we write equation 
5. 36 in dimensionless form, replacing the finite "delta" quantities 
with differentials^ then substituting 3. 15 and 5. 35, we have 
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cl i - 



ITr 

Vt. 



di 



P 0 /z U-e) 



?/z 



(6.9) 



where again, di is in effect one. Substituting equation 6. 4 for 
d 1 and. neglecting "end effects" in the integration, we have 



4* s ITe* 



-/ 



)iz - ^oe-ize 2 -] 



n ■ ,£9si*e -+■ , , 

’ 3 60 r fl-e* 



(6. 10) 






The "end effects" we mention are the initial cutoff angle and the final 
cutoff angle not being equal so that we have a small amount of time 
neglected for each cycle due to the shift of the cutoff angle; the final 
energy is added before starting the final coasting orbit; and the initial 
energy is added during the latter part of the initial coasting orbit* 

This might be clearer if we note that the integration gives us 



T - 4?. + A L 



*<■ 



-h- 



+■ aT. 



’*+/ 



-■f-i 



+ A 7 



e-f 



( 6 . 11 ) 



where /l "Tgy ^e trave ^ ti me ^ or 2 // radians at eccentricity 

Qj , whereas the true mission time would be given by: 



T-AT.^AT +dT + aT +JT-+'-' + dT +AL + * t7 x,. 

1 <?/,, t+t tv* * f V-/ f 1 



( 6 . 12 ) 



where is the thrust time for the initial cycle, 

o . — ~ n & . , and dT- is the cutoff angle shift 

e ] ~ e f., J 

time. If the number of cycles is not large, then equation 6. 10 should 
be corrected for "end effects" by adding 

e z + --k 



T = AT. I l -h -&■ o 

CoKR i l Q -f- fe 



-at: - AT. 






f 



( 6 . 13 ) 
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where is the initial cutoff angle given by equation 6 . 2 » Q-f 

obtained by evaluating equation 6. 3 at <2 = £ f . The cutoff angle 
shift time correction is an average value obtained by using the first 
term of equation 6. 1. The other corrections follow simply by com- 
paring equations 6. 11 and 6. 12. 
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VII. RESTRICTED PERIGEE CYCLE WITH 



BATTERY STORAGE 



We must still consider the possibility of different values of 
specific power plant mass, K > for the two systems, since we have 



the capability of storing energy with the discontinuous method' while 
coasting if we provide batteries* The energy generated by the power 
supply while coasting could be stored in batteries and then used in 
addition to the power supply while thrusting. 



pack is included with the power supply and that other conditions of 
the thrust program are unchanged; later the battery will be consid- 
ered as part of the payload* We shall further assume that during the 
early part of the program, when £ is small, that the coasting time 
is not sufficient to fully charge the batteries, but as £ increases, 
a point is reached where the batteries do become fully charged, and 
are fully recharged during each subsequent coasting period. We 
have equation 2. 18, evaluated at the “burnout condition 11 : 



consider P as a function of time* When the batteries are not fully 
charged, the power is given by 



To study this possibility, we shall first assume that a battery 




(7. 1) 



where we 



have multiplied and divided by (K/M^) , and must now 



p — A Tc vc l 



K AT 



(7. 2) 



neglecting battery losses while charging, where is the mass of 
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VII. RESTRICTED PERIGEE CYCLE WITH 



BATTERY STORAGE 



We must still consider the possibility of different values of 
specific power plant mass, K > for the two systems, since we have 



the capability of storing energy with the discontinuous method: while 
coasting if we provide batteries 0 The energy generated by the power 
supply while coasting could be stored in batteries and then used in 
addition to the power supply while thrusting* 



pack is included with the power supply and that other conditions of 
the thrust program are unchanged; later the battery will be consid- 
ered as part of the payload* We shall further assume that during the 
early part of the program, when €. is small, that the coasting time 
is not sufficient to fully charge the batteries, but as C increases, 
a point is reached where the batteries do become fully charged, and 
are fully recharged during each subsequent coasting period. We 
have equation 2. 18, evaluated at the "burnout condition' 1 : 



To study this possibility, we shall first assume that a battery 




( 7 . 1 ) 



where we have multiplied and divided by and must now 

consider P as a function of time* When the batteries are not fully 



charged, the power is given by 



p — ^ Tc ycue 



K AT 



(7. 2) 



neglecting battery losses while charging, v/here 



is the mass of 
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the electrical generating unit, K is the specific mass of the gener- 
ating unit, AT is the constant thrusting time per cycle, and 
A T C y C ^ e is the travel time for S7T radians given by equation 
5„ 3 6. The batteries will be fully charged when 

< 7 ' 3; 

" batt. ^ 

or 



^orbit^ Ar+ 



M 



batt. 



K 






K 



(7.4) 



J lx batt. 

are the mass and the specific mass of 
the battery, respectively. From equation 7. 1, we define , cor- 

responding to 2. 22: 



where M, ,, and K, 

batt. batt. 



0 



-if T H 

2- 1 P 

6 



(7. 5) 



But since £ is constant, and 



Y‘ = 



_ \ 



r 



Mad + 



I- 



- Ml^ -h Mda.fi-. : 

) c/t ) 

a* fil’ 



(7.6) 



where -f (£) is introduced to remind us that the integrand is 
zero when not thrusting. Using equation 7. 4 in 7. 2 , we can write 
the maximum power available when the batteries are fully charged: 




K l 

Mf !<ui+. XT 



(7.7) 
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Using equations 7. Z and 7. 7 in 7. 6 : 






iC/- IT. 

■fit) &T 

L'o ^ 1 °™^T 



ol-t -f 




(7. 8) 



where *t . is the time at which the batteries become fully 
c r i Lo 

charged* Now, the integrals in equation 7* S may be treated as sum- 
mations since the thrust is on for the constant time A'T' during 
each cycle: 






r 



2 _ 



£1 

z 



K(l 



+ 




(N-M cieir ;zir 
+ J_ j< \ 

KuJ J 



(7. 9) 



where f\l .. is the cycle number after which the batteries are 
'mrit. 

fully charged. Using equations 6. 9 and 7. 4, we can solve for the ec- 
centricity when the batteries become fully charged: 



F? 



Zir 



AT + 



NW K 



Then, using equations 6. 4 and 6. 9 : 



Ncnr No- it 

Z ^ II 6 r 

* H 




(7. 10) 
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3££ ( i^e) 3/a 




(7. 11) 



Integrating this equation and determining N - IV crit 



from equation 
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6. 5, using 7. 10, equation 7. 9 becomes: 

1? P/f. 



2 ^ 2 R’ 4 



£ l 9 l 

y 2 = IJi )</, 

J 2 JM 4€fe( V 277 



(6T> f (2- |)F?-2) 
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7ga/ f l8Z8e-57 / >e Z + ZlC€? p* 1 



/, , x,i , 38Zl + l3Z8e-5 ??< 
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J e, 



(7. 12) 



where 



I + P 

and 

| -h «* 



^ baft. 



_J f2 b ath K 

AT M ? |<utt. 



(7. 13) 



(7. 14) 



Now if we form the ratio of equation 2. 22 for the spiral orbit and 
equation 7. 12, multiplied by the factor I 4 p , we have 
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T Vf 1 (€j/£ s ) P 









2 77 



-Ke- r)lfi r ?'-2 : ) + 



Qg£f-230et -h\063e c -2404 gf -3</<e c + 4gJ4 

e ll i -ef 



f i~SlA/e c - g-FP.7 ^ -fe* 

F 3 J? 6 0 




-f 



4 - ?llLLl^£z£21^it ^ lge? 
1 T/TTe^ 




(7. 15) 



When F"Ry 2 is less than one, the spiral method of energy addi- 
tion must have a Y less than the Y for the restricted perigee 
method, and by equation 2. 25 the spiral method must have a greater 
payload. If FRyZ is greater than one, then there may be a pos- 
sibility of carrying a greater payload by using the restricted perigee 
method. This will depend on the ratio of K/Kfoatt that can be 



obtained. 
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VIII. COMPARISON OF RESULTS 

Now to compare the results we have obtained, equation 5. 31, 
the spiral thrust time, is plotted against the energy added, equation 
5. 32, in fig. 3. Using a value of ft - 0. 8 , which corresponds to 
a maximum ©^ = 0. 944 by equation 6. 3, the restricted perigee 
thrust time is computed using equations 6. 2, 6. 5, and 6. 6 for values 
of eccentricity from zero to 0. 95, with = 0. 95 . The corre- 

sponding increase of energy is computed using equation 6. 7, and 
these results are plotted in fig, 3, Since the total time must be equal 
for either method at a specified energy addition, we may compute 
equations 5. 31 and 6. 10 at specific energy levels and take their ratio 
to plot /£$) in fig* 3, This is the ratio of acceleration for 

the discontinuous to the spiral that is necessary to accomplish the 
mission. 

It seems clear from fig. 3, in view of the acceleration ratios 
required for the mission and equations 2. 13, 2. 22, and 6. 8, that T 
will be greater for the discontinuous method if K is the same for 
both methods and starting from the circular orbit with ” 0. 95 , 

For these conditions, the discontinuous thrust method will provide 
less payload. Further, if we decrease A , the initial cutoff angle, 
we see from equation 6. 10 that £ must be increased for the same 
total time; thus, Y will be even larger for the discontinuous 
method, resulting in even less payload. 

Using equation 7. 10 and fig. 3, we have evaluated equation 
7. 15 for three different final energies, expressed as eccentricity, for 
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various values of o< . These results are plotted in fig. 4. 

From reference (3), we find that a typical value of K would 
be about 150 pounds per kilowatt. This value will decrease as space 
power plant technology progresses. In reference (3), we also find 
^ batt. of 75 pounds per kilowatt-hour for nickel-cadmium batteries, 
which have very reliable charge and discharge characteristics, and 
14 pounds per kilowatt-hour for silver-zinc batteries, which are not 
as reliable for a large number of recharging cycles. From equation 
7. 3, we see that both values of Kh must be divided by 

O cL L C© 

Jr 

3600 (^ 0 / R 0 ) to have the proper non-dimensional time units. 

Thus, we have the values for of 16. 3 for nickel -cadmium 

and 3. 13 for silver-zinc. 

We cannot improve over the spiral thrust program with either 
battery supply if the final energy must correspond to an eccentricity 
of 0. 95, since in fig* 4 the maximum PRy2 is 1. 006 . 

Using equations 5. 35 and 6. 9 3 we compute the non-dimension- 
al thrust time for one cycle, corresponding to an initial cutoff angle 



of 0. 8 and 



€ 0 = 0 



* 7 



_ (Z)(0-B) d - _ I- _ j ’-j g Q 



(o.95) Vi 



( 8 . 1 ) 



Z1T 

Using this and equation 7. 14, we compute the values of JA for the 
nickel-cadmium battery, for <=* = 7. 5 and c< = 10. 0, respectively: 






= 1.4-S 






= 1.33 



( 8 . 2 ) 



and for the silver- zinc battery: 
fA = 0.2 7 , ^ - 0-3/0 



(8.3) 
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From fig. 4, we read FRya^l. 21 for the final eccentricity of 

0. 4 and cx = 7. 5 . But from equations 7. 15, 8. 2, and 8. 3, the 
values of JA for either battery will require oj to be greater 
than Y s ) again, we find no advantage in the restricted perigee 
method of energy addition. However, for the final eccentricity of 0. 7 
we read from fig. 4 FRyZ of 1. 29 and 1. 39 for o< *s of 7. 5 and 
10. 0, respectively. From equations 7. 15 and 8. 3, we compute for 
this condition: 



- 1.00 8 , -£-= 1.010 < 3 - 4 > 

Yj > Yj 

for the two values of o< , respectively. From equation 2. 25, we can 
see that in this case the restricted perigee method has a slightly 
larger payload than the spiral method for this particular final energy. 
To see this result more clearly, we form the ratio of equation 2. 25 
for the two methods using 8. 4: 

(M ‘ + = i 1 - lit (8.5) 

(M, 4 Ms) s (I -l.onjf 

for o< = 10. 0 . Now, clearly, the numerator of 8. 5 is greater than 
the denominator, so the ratio must be greater than one. 

Although we have found that the restricted perigee method of 
energy addition, using silver-zinc batteries, yields a slightly larger 
payload than the spiral method for certain values of final energy, we 
still have the problem of battery reliability. There is no question of 
the recharging reliability of the nickel -cadmium battery, but we 
found the specific mass, K - Q&tt > to he too high to be useful. The 
silver-zinc battery could be used for a slight payload advantage, be- 
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cause of its lower specific mass, but its recharging reliability must 
first be improved. 

We must finally compare the two methods of energy addition 
when reliable, rechargeable batteries are part of the payload but are 
available to be used in the thrust program. In this case, JA = O and 
we see immediately from fig. 4 that there are certain values of 
and final energy where the restricted perigee has a decided advantage 
over the spiral method. 

As an example, let us determine the advantage where a 
nickel-cadmium battery, K-^ ^ ~ 16. 8 , is part of the payload. 

Using equations 7. 14 and 8. 1, we compute values for 
for cx = 7. 5 and o< = 10. 0, respectively: 

Fl batt. , M batt 



Me 



1. 45 



Ml 



= 1. 93 



( 8 . 6 ) 






Even though the ratio of battery to power supply mass is rather 
large, it is not unrealistic since we could easily have a satellite 
where both battery and power supply are but a small fraction of the 
initial mass. From fig. 4 we read FR = 1° 21 for the final 

eccentricity of 0. 4 and = 7. 5 , and since ju = O : 

= 1.10 < s - 7 > 
h 

Forming the ratio of equation 2. 25 for the two methods, using 8. 7 : 



_ Cl -Yj f 

(M, t-MsJj (l-IAOYjf 



Similarly, from fig. 4 we read FR yz of 1. 29 and 1. 39 for the 
final eccentricity of 0.7 and o< 's of 7. 5 and 10. 0, respectively. 
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Since - O , we have immediately: 

-~ = I.13G , -^= 1.179 (8.9) 

•d x d 

for the two values of <=< respectively. Again forming the ratio of 
equation 2. 25 for the two methods, using 8. 9> we have 

_ (I ~^) 2 (m, + m s )j (i-rj) 2 

(M, + Ms)s (I-I.13G Yj) z ) (M, (J -M?3^() a (S * 10) 

for the two values of = 7. 5 and = 10. 0, respectively* The 

ratios in equations 8. 8 and 8* 10 will depend on the value determined 
for > computed from equation 7. 12, but they do show a definite 

advantage for the restricted perigee method of energy addition. 

For example, suppose the value for is 0* 25, that is, 

let 25 per cent of JM C be propellant* Substituting in equation 8* 10, 
we compute ratios of 1* 097 and 1* 131 for the two values of o< , re- 
spectively. Now, this is about a 10 per cent gain in payload and 
structure mass by using the discontinuous rather than the continuous 
method of thrust programming, where the batteries are part of the 
payload* 

It seems rather surprising that v/e have found no advantage 
for the restricted perigee method of energy addition for the final ec- 
centricity of 0* 95, where we found a definite advantage for the final 
eccentricity of 0. 7. There are two reasons for this result* First, 
we can see from equation 6* 9 that the travel time for 2.7 T radians 
increases as the eccentricity increases. Thus, if we thrust for the 
same amount of time during each cycle, there is a greater percentage 
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of coasting time as the eccentricity increases,, This requires an in- 
crease of (^-d/^s) 9 as s h° wn i- n fig. 3, to accomplish the 

mission* and in turn, we find a decrease in the ratio of (Y s /Yd) 
so that there is no advantage when reaches 0* 95 . Second, by 

restricting the cutoff angle to one radian or less, we do not take full 
advantage of the greater velocities provided by the elliptic orbit as 
the eccentricity approaches one* This can be seen clearly in fig. 2 
where, at a cutoff angle of one radian and <8 = 0* 95 , the specific 
power of the restricted perigee is many times greater than the spiral 
method of energy addition* To obtain the full advantage of this much 
greater specific power, we must go to larger cutoff angles* 

A study of the restricted perigee method of energy addition 
at larger cutoff angles would require the use of numerical methods. 
This would be rather involved, using equation 4* 13 to evaluate A <£ 
for various values of Q 7 and & , then using this result for a 

corresponding evaluation of the change in energy, equation 4. 12. 

After tabulating these "cycle 11 results, summations of time and energy 
could be made with a corresponding calculation of for various 

values of oc . By making such a numerical study, we can clearly 
define the area where the restricted perigee method, with batteries, 
has an advantage over the spiral method of energy addition* We will 
also be able to determine if there is any area where an advantage ex- 
ists without the use of batteries* 
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IX. CONCLUSION 

We have found that the discontinuous method of energy addi- 
tion, using a restricted perigee, has no payload and structure mass 
ratio advantage over the spiral method when a specified energy is 
imparted in a specified time and the cutoff angle is restricted to one 
radian or less. There is the possibility of a very small advantage, 
for a certain range of final energies, if silver-zinc batteries are in- 
cluded in the power supply to store energy during the coasting periods, 
providing their reliability can be improved to an acceptable level. On 
the other hand, a moderate decrease of the power-plant specific mass 
will rule out any advantage of including batteries in the power supply. 

In the case where batteries are included in the payload, say, 
to periodically energize a transmitter after the satellite has been 
placed in final orbit, we can definitely increase the payload and 
structure mass ratio for certain values of final energy by using these 
batteries in the restricted perigee thrust program. This does not ap- 
ply in general, obviously, since the mass fraction of the batteries re- 
quired in the payload must be compatible with the mass fraction of 
batteries required for the thrust program to accomplish the mission.. 
When we do have a compatible battery requirement, then we can in- 
crease the payload and structure mass ratio by using the restricted 
perigee method of energy addition. 

Since we are restricting our investigation to low thrust de- 
vices, a numerical study of the restricted perigee method of energy 
addition at cutoff angles greater than one radian could be made using 
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the perturbation equations developed earlier. Such a study would 
follow the procedure briefly outlined in the last paragraph of the pre- 
ceding section, and is necessary to clearly define the area where the 
restricted perigee method can be used to advantage. 
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